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static  Cavity  Expansion  Model  for  Partially  Confined  Targets 


Yehuda  Partom 


Abstract 

The  cavity  expansion  model  (CEM),  originally  proposed  as  an  indentation  theory,  has  been 
used  extensively  to  estimate  the  resistance  of  targets  to  long  rod  penetration.  In  the  classical 
model  the  target  is  infinite  and  the  cavity  is  opened  up  from  zero  radius.  In  previous  work  we 
applied  the  classical  approach  to  laterally  finite  targets  and  we  ran  into  difficulty.  The  resis¬ 
tance  came  out  as  decreasing  with  cavity  radius.  To  circumvent  the  difficulty  we  introduced 
an  averaging  scheme  with  a  free  parameter  to  be  determined  from  computer  simulations. 

In  this  work  we  use  a  different  approach,  we  open  up  the  cavity  from  a  Hnite  radius.  We  first 
propose  the  concept  that  the  target  penetration  resistance  is  the  limit  of  stability  of  the  cavity. 
We  then  apply  this  concept  to  estimate  the  resistance  of  partially  confined  targets.  We  do  this 
for  both  cylindrical  and  spherical  cavities. 

Comparing  the  model  solutions  to  computer  simulation  results  reported  in  [2]  we  find  good 
agreement  for  the  spherical  CEM,  but  not  for  the  cylindrical  CEM. 


1.0  Introduction 

In  [1]  we  applied  the  static  cavity  expansion  model  (CEM)  in  cylindrical  symmetry  to  estimate 
the  efficiency  of  lateral  self-confinement  in  metal  and  ceramic  targets.  In  [2]  we  calibrated  and 
validated  the  results  obtained  in  [1]  using  computer  simulations. 

Our  approach  in  [1]  and  [2]  was  along  the  lines  of  the  original  CEM  [3]  namely,  evaluating  the 
internal  pressure  P(a)  needed  to  open  a  cavity  from  zero  radius.  For  a  finite  outside  radius,  b, 
we  found  that  P(a)  was  not  constant  but  decreased  with  the  cavity  radius  “a.”  We  therefore 
averaged  P(a)  over  a  certain  range  of  “a”  and  defined  this  average  to  be  the  target  resistance. 
We  calibrated  the  averaging  range  from  the  results  of  computer  simulations. 

In  what  follows  we  use  a  different  approach.  We  first  show  that  the  value  of  P(a)  obtained  from 
the  classical  CEM  (b  >  <»)  is  actually  the  limit  of  stability  of  a  finite  radius  cavity  expanded 
by  internal  pressure.  We  then  use  this  insight  and  evaluate  the  limit  of  stability  for  finite  b.  We 
regard  the  limit  of  stability  P^  to  be  the  target  penetration  resistance. 
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We  find  that  depends  on  bo/oQ  (where  a^,  bQ  are  the  initial  inside  and  outside  radii).  Identi¬ 
fying  ao  with  the  projectile  radius  we  obtain  good  agreement  with  computer  simulations  for  a 
spherical  cavity  but  only  qualitative  agreement  for  a  cylindrical  cavity. 

In  Section  2  we  explain  the  concept  of  the  limit  of  stability.  Section  3  provides  the  analysis  for 
a  cylindrical  cavity,  and  Section  4  provides  the  analysis  for  a  spherical  cavity.  Chapter  5 
shows  the  model  solutions  for  penetration  resistance  as  a  function  of  the  degree  of  confinement 
and  compares  this  to  computer  simulation  results. 


2.0  Limit  of  Stability 

To  explain  the  limit  of  stability  concept  we  evaluate  the  elasto-plastic  cavity  expansion  from  a 
finite  initial  cavity  Uq. 


2.1.  Cylindrical  Cavity 

The  elastic  field  is: 


pY  c 
u  =  ^  — 
2G  r 


Gr-Oe 


where: 

u  =  radial  displacement 

CT^,  Gq  =  radial  and  tangential  stress  components 
P=  V3/3 

Y  =  yield  stress  and  flow  stress 
G  =  shear  modulus 
c  =  elastic-plastic  boundary  radius 
r  =  radial  coordinate 


(1) 
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and  where  at  r  =  c  we  have: 


a  -a0  =  -2pY, 

which  follows  from  the  von-Mises  yield  surface  and  from  [3]: 
Gz  =  ^(a^  +  ae)  , 

where  a  =  axial  stress  component.  At  r  =  c  we  have; 


^  ^  PY 

“(c)  =  fee  , 


a,(c)  =  -PY. 

From  mass  conservation  (ignoring  density  changes)  we  have: 
c^  -  a^  =  [c  -  u(c)]^  -  ao^, 


a^'PY 


The  plastic  field  is: 


V 


o,  =  2pY(-i.l„l), 


and  the  cavity  wall  pressure  is: 


P(a)  =  PYh+ln 


pY 


f  a 

V  a  yj 


We  see  that  as  ao  goes  to  zero,  P(a)  increases  to  a  limit  value: 


Pc=  P(a)  =pY(l+ta-^l 

PY 

ao^O  ' 


which  is  the  classical  result  for  cylindrical  cavity  expansion. 


The  P(a)  curve  from  (7)  is  shown  in  Figure  1  and  discussed  below  together  with  the  spherical 
cavity  results. 


2.2.  Spherical  Cavity 

The  elastic  field  is: 


u 


Y 


3 

C 


6G  J.2’ 


3 


r 


where  (r,  6,  (p)  are  the  spherical  coordinates. 
At  r  =  c  we  have: 


u(c) 


a,(c)  =  -?Y, 

(O.-Gq)  (c)  =  -Y. 

From  mass  conservation  (and  no  density  change)  we  have: 
c^  -  a^  =  [c  -  u(c)]^  -  ao^, 


2G 

Y 


(9) 


(10) 
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The  plastic  field  is: 


(12) 


and  the  cavity  wall  pressure  is: 


P(a)  =  1+ln 


Y 


(  „  3\-| 


1- 


^0 


V  ay 

from  which  we  get  the  limit  value: 


(13) 


P,  =  P(a)  =|Y('l+ln^]. 

(ao  0) 

The  P(a)  curve  from  (13)  is  also  shown  in  Figure  1. 


(14) 


2.3.  Limit  of  Stability  and  Penetration  Resistance 

The  cavity  expansion  results  obtained  for  QPi  =  100  are  shown  in  Figure  1. 


Figure  1.  Cavity  expansion  (a/ao),  in  an  inHnite  domain,  as  a  function  of  scaled  internal 
pressure  (P(a)/Y),  for  cylindrical  and  spherical  symmetry.  The  horizontal  asymptotes  are 
the  limits  of  stabUity  of  the  cavity.  is  given  by  (7),  P^  by  (13),  P.,  by  (8)  and  P«,  by  (14). 
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We  see  that  as  long  as  P(a/ao)  <  P^,  the  cavity  is  statically  stable  with  a  radius  “a.”  At  P  =  Pj 
the  cavity  can  grow  indefinitely,  quasistatically  (with  zero  velocity).  Also,  because  a/aQ  is  in¬ 
finite  when  a  finite  cavity  is  opened  from  zero  radius,  P5  is  identical  to  the  penetration  resistance 
obtained  from  the  classical  CEM  [3].  This  provides  a  different  way  of  looking  at  penetration 
resistance  obtained  from  CEMs:  penetration  resistance  is  the  internal  pressure  on  the  cavity 
wall  that  would  make  the  cavity  expand  indefinitely,  at  zero  velocity.  For  P  >  P5,  the  cavity 
wall  is  accelerated  by  the  internal  pressure. 

In  [1]  and  [2]  we  used  the  CEM  to  estimate  the  penetration  resistance  of  partially  confined  tar¬ 
gets.  Using  the  classical  approach  (opening  a  finite  cavity  from  zero  radius)  we  found  that  the 
internal  pressure  needed  to  open  the  cavity  to  a  radius  “a”  was  a  decreasing  function  of  “a.”  To 
extract  a  penetration  resistance  value  from  P(a)  we  averaged  it  over  the  range  0  <  a  <  a^  where 
aj.  was  calibrated  from  computer  simulation  results. 

By  applying  the  limit  of  stability  approach  to  cavity  expansion  in  a  finite  domain  we  can  avoid 
the  arbitrariness  of  the  averaging  approach. 


3.0  Cylindrical  Cavity  Expansion,  Finite  aQ  and  bo 

The  elastic  field  is: 


u 


pYfc^  G 

20**  2'^X  +  g 
Vr 


a,  =  -PY 


C  C 
I  ~  ”^2  I 

^r  b  J 


ae  =  PY 


o.  =  PY 


c 


\*a  ^2- 


=  -2PY 


P  =  -PY 


K  c 


A,  +  G 


(15) 
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where: 


X  =  Lame’s  modulus, 

K  =  Bulk  modulus. 

The  elastic  moduli  ratios  in  (15)  are: 


_L_ 

A,  +  G 


=  2v, 


G 

A,  +  G 


l-2v. 


K 

X  +  G 


|(1+V), 


where  v  =  Poisson’s  ratio, 
and  at  r  =  c  we  have: 


^  ^  PY 
u(c)  = 


l  +  (l-2v)^ 


a,(c)  =  -pY 


The  plastic  field  solution  is: 


= 


-PY 


^  2  2^ 

1  C  ,  C 

1  — In— 

V  2  2 


V 


r  J 


so  that: 


P(a)  =  pY 


f  2  2\ 

1  C  .  c  ' 

1  — 7  +  In— 

I  u2  2 

V  b  ay 


From  mass  conservation  we  have  (as  in  Section  2.1): 


(16) 


(17) 


(18) 


(19) 
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(20) 


=  [c  -  u(c)]^  - 


a  J 


c 

u(c)’ 


G 

PY 


2  2 

l-aoVa 


l  +  (l-2v)-^ 
b 


(21) 


Also,  at  r  =  b  we  have: 
b  =  bo  +  u(b). 


u(b)  =  4&'-v). 


so  that: 


c 

a 


^  ^  c 
ao  a  b 


•(1-v)-^ 

b^ 


(22) 


(23) 


Equations  (21)  and  (23)  are  nonlinear  simultaneous  equations  in  c/a  and  c/b.  Substituting  c/a 
from  (23)  into  (21)  we  obtain  a  quadratic  equation  in  c^/b^  which  can  be  solved  analytically  for 
given  bo/ao  and  a/aQ.  Finally,  substituting  into  (19)  gives  P(a)  for  given  bo/ao  and  a/aQ. 

When  the  elastic-plastic  interface  reaches  the  outside  boundary  we  have: 

c  =  b. 


P(a)  =  PYln^, 


(24) 


b^  _  G  1-aoVa^ 
^2  -  Py’  2(1 -V)  • 

But  from  (22): 


8 


b 

b 


1 


(25) 


Eliminating  b  from  the  last  of  (24)  and  (25)  we  get: 


-  =  1+  - 


b.A2  2(1 -V) 


G 


'l-(l-v) 


G. 


for  c  =  b. 


(26) 


For  a  >  ac=b  the  second  of  (24)  still  holds,  and  b^/a^  is  obtained  from  mass  conservation: 


b^  -  a^  =  bo^  -  ao^, 


^  =  1+-? 
2  2 
a  a 


2/u  2  A 


(27) 


Vao  ) 


Results  for  P(a)Ar  as  a  function  of  a/aQ,  for  GAf  =  100,  v  =  0.3,  and  for  different  values  of  bo/aQ, 
are  shown  in  Figure  2.  We  see  that  the  curves  rise  to  a  peak  and  then  drop.  The  value  of  P(a) 
at  the  peak  is  the  limit  of  stability  pressure  P{.  We  see  that  Pf  increases  with  b(/ao  and  that  the 
corresponding  a/ao  value  increases.  In  Chapter  5  we  further  discuss  these  results  together  with 
the  spherical  cavity  results  obtained  in  Chapter  4. 
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0 

I  1  1.5  2  2.5  3  3.5  4  4.5  5 

I  a/ao 

Figure  2.  Cylindrical  cavity  expansion  (a/ao),  ^  domain,  as  a  function  of  the  scaled 

internal  pressure  (P(a)/Y),  for  G/Y  =  100  and  v  =  0.3.  The  different  curves  are  for  different 
values  of  bo/ao* 
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4.0  Spherical  Cavity  Expansion  (finite  Oq  and  bo) 


The  elastic  field  is: 


_  Y  ,  o  1  -  2v 

“  -  7  +  ^  ■ 

vr  D  y 


^  3  3A 

a  =  --Y  — -  — 

>■  3  3  3  ’ 

[r  b  J 


_  _  2  f  1  ^ 

®  n2‘,3  ^^3!’ 


a,-ae  =  -Y-^, 


p  =  -2y.£: 
3  b^’ 


and  at  r  =  c  we  have: 

u(c)  =  ^cfl  +  2-i2:^.^ 

6G  1  +  V  1^3 


a,(c)  =  -?Yl-L 
V  b 


The  plastic  field  solution  is: 


r,  f  3  3\ 

ctr  -  --Y  1  -  —  +  ln—  , 
V  b  r 


so  that: 


^  f  3  3 

P(a)  =  3-Y 

V  D  a 
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From  mass  conservation  we  have: 


-  a^  =  [c  -  u(c)]^  -  ao^ 


^  a^ 
V  ay 


c 

u(c)’ 


3 

C_ 

3 


a 


2G 

Y 


(1-aoV) 


1+2- 


l-2v 
1  +  V 


Also,  at  r  =  b  we  have: 
b  =  bo  +  u(b). 


Y  c 

u(b)  = 

^  ^  2G  ^^2 


1-v 
1  +v’ 


so  that: 

^  ^  c 
c  _  ap  a  b 

^  ,  Y  1-v 

2G' 1+v’^ 


(32) 


(33) 


(34) 


(35) 


Equations  (33)  and  (35)  are  nonlinear  simultaneous  equations  in  c/a  and  c/b.  Substituting  c/a 
from  (35)  into  (33)  we  get  a  cubic  equation  in  c^^^  which  can  be  solved  numerically  with  a 
standard  bisection  routine.  Finally,  substituting  back  into  (31)  gives  P(a)  for  given  bp/aQ  and 
a/ao- 

When  the  elastic-plastic  interface  reaches  the  outside  boundary  we  have: 
c  =  b, 

3  1  3/3 

b-^  _  l-ap/a 

^  "  3  Y  1-v’ 

2  G'  l-HV 

9 

P(a)  =  ^Yln^. 

a 


11 


But  from  (34): 


b  _  1 

bo  “  i.X.izv' 

2G  1+v 

Eliminating  b  from  the  second  of  (36)  and  (37)  we  get: 


1+1- 
^0 


3  Y  1-v 
2G’  1+v 


1- 


_Y  1-v 
2G  ■  1  +  V 


for  c  =  b. 


(37) 


(38) 


For  a  >  a^-fj  the  third  of  (36)  still  holds,  and  b^/a^  is  obtained  from  mass  conservation: 
b3-a3  =  bo3-ao3, 


u3  /a  ^ 


3  ^ 


(39) 


Vao 


y 


Results  for  P(a)nr  as  a  function  of  a/a^,  for  G/Y  =  100,  v  =  0.3,  and  for  different  values  of  bo/ao, 

are  shown  in  Figure  3.  The  curves  are  similar  in  nature  to  those  shown  in  Figure  2  for  the  cy¬ 
lindrical  cavity. 


i  1  1.5  2  2.5  3  3.5  4  4.5  5 

i  a/ao 

Figure  3.  Spherical  cavity  expansion  (a/a^,),  in  a  finite  domain,  as  a  function  of  the  scaled 
internal  pressure  (PfalA^),  for  G/Y  =  100,  v  =  0.3.  The  different  curves  are  for  different 
values  of  bo/ag. 
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5.0  Penetration  Resistance  Results  and  Conclusions 

As  explained  in  the  introduction,  the  peaks  of  the  curves  in  Figures  2  and  3  are  the  limits  of 
stability  of  the  cavity,  in  terms  of  the  scaled  internal  pressure  (P^/Y).  In  Figure  4  we  show 
as  a  function  of  bo/a^  for  cylindrical  and  spherical  cavities. 


Figure  4.  Limit  of  stabUity  (maxima  of  P(a)/Y  versus  a/a^  curves)  as  a  function  of  bo/aQ 
for  cylindrical  and  spherical  cavities.  G/Y  =  100,  v  =  0.3.  Asymptotes  are  the  same  as 
in  Figure  1. 


We  see  that  the  spherical  cavity  curve  approaches  its  asymptote  faster  that  the  cylindrical  cavity 
curve.  In  Figure  5  we  show  that  the  same  results  only  normalized  to  their  asymptotic  values 
(Pf  for  bo  oo).  We  also  identify  P^  as  the  penetration  resistance  R  and  denote  the  asymptotic 
values  by  Rjjjf. 


normalized  penetration  resistance 


Figure  5.  Same  as  in  Figure  4  but  normalized  to  values  at  infinite  also,  is  identified 
with  penetration  resistance  R.  Points  are  from  simulations  with  steel  targets  in  [2]  as¬ 
suming  that  D/Dp  =  bo/aQ. 
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We  also  put  on  Figure  5  the  results  of  numerical  simulation  with  steel  targets  from  [2],  In  doing 
that  we  assumed  that  D^/Dp  in  the  simulation  (where  Dj  =  target  diameter  and  Dp  =  projectile 
diameter)  is  identical  with  bo/aQ  in  the  CEM.  We  see  that  the  simulation  points  are  quite  close 
to  the  spherical  cavity  curve.  (For  steel,  G/Y  =  100  and  v  =  0.3  is  a  good  approximation,  and 
the  sensitivity  to  those  parameters  is  low.) 

We  conclude  that: 

•  The  limit  of  stability  approach  is  an  appropriate  tool  for  estimating  penetration  resis¬ 
tance  of  partially  confined  targets. 

•  The  spherical  CEM  is  better  suited  to  model  penetration  resistance  than  the  cylindrical 
CEM. 
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